6-84  ʹnlpredciģԤ⡣
>>  load reaction
>> [beta,resids,J] = nlinfit(reactants,rate,@hougen,beta);
>> [ypred,delta] = nlpredci(@hougen,[100 300 80],beta,resids,J)
ypred =
   10.5940
delta =
    0.4456
>>
